Routes to thermodynamic limit on scale-free networks 
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We show that there are two classes of finite size effects for dynamic models taking place on a 
scale-free topology. Some models in finite networks show a behavior that depends only on the 
system size N. Others present an additional distinct dependence on the upper cutoff k c of the 
degree distribution. Since the infinite network limit can be obtained by allowing k c to diverge 
with the system size in an arbitrary way, this result implies that there are different routes to 
the thermodynamic limit in scale-free networks. The contact process (in its mean-field version) 
belongs to this second class and thus our results clarify the recent discrepancy between theory and 
simulations with different scaling of k c reported in the literature. 
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The thermodynamic limit is a crucial concept in sta- 
tistical physics. For example, no true singularity can 
occur in a system composed of a finite number N of ele- 
ments at nonzero temperature so that the very con- 
cept of phase-transition is well defined only in the limit 
A" — > co. However, all real systems are finite, and more 
so in numerical simulations, and to study their proper- 
ties one has to understand the role of finite size effects 
when the limit N — > co is taken. In the case of critical 
phenomena, the theory of finite size scaling (FSS) [2| has 
successfully accomplished this task for processes on reg- 
ular lattices, allowing to detect the signature of continu- 
ous phase-transitions even in very small systems. Finite 
size effects are expected to be all the more relevant for 
systems with a strongly heterogeneous interaction pat- 
tern, such as those described in terms of complex net- 
works 0, Hj]. Indeed, because of the small- world prop- 
erty [H observed in most networks, the number of neigh- 
bors that can be reached starting from a certain node 
grows exponentially or faster with the geodesic distance. 
This implies that, even for large networks, just a few 
steps are sufficient to probe the nniteness of the system. 
Moreover, most real networks exhibit the scale-free (SF) 
property i.e. they have a degree distribution P(k), 
defined as the probability that a vertex is connected to 
other k vertices, decaying for large k as P(k) ~ fc -7 with 
2 < 7 < 3, so that the local topological properties dis- 
play very strong fluctuations, increasing with the size of 



the network 12j. In order to understand phase transi- 



tions, and in general any kind of dynamical process, in 
SF networks it is thus necessary to first comprehend how 
size effects, and in particular FSS, work in this class of 
systems. 

Some attempts in this direction have been already 
done. For equilibrium continuous phase transitions the 
situations seems to be rather well established in some 
cases. For example, a FSS phenomenological theory 
has been developed for the Ising model [1, 0], based 
on a scaling ansatz for the free energy, leading to the 



FSS form for the magnetization at zero external field 
m(A,N) = Af-^/fAJV 1 /"), where A is the distance 
to the critical point, and the critical exponents (3 and v 
depend on the degree exponent 7 and are defined only for 
7 > 3 22]. In the context of nonequilibrium phase tran- 
sitions, however, the situation is not so clear. In the case 
of the contact process (CP) in SF networks Q Hong et al. 
[H proposed a FSS for the average density of particles in 
surviving runs scaling as p s (A,N) = N-?/" /(AN 1 ^), 
with (3 = 1/(7 - 2), v = (7 - l)/(7 - 2) for 2 < 7 < 3, 
and j3 = 1, v = 2 for 7 > 3. This theory, however, is 
still under debate [9|, [10( , since numerical simulations in 
random neighbors SF networks yield incom pat ible results 
0. For other kinds of FSS, see also Ref. [111. 

In this Letter we provide a step forward in the under- 
standing of finite size behavior in dynamical processes 
on SF networks by showing that two different scaling 
scenarios may occur. A first class of processes exhibits 
finite size effects depending exclusively on the network 
size N. On the other hand, a second class of processes 
displays anomalous finite size effects, in the sense that 
their properties depend explicitly and independently not 
only on the number of nodes N, but also on the upper 
cutoff k c of the degree distribution [12|. Hints to this 
fact can be found in previous works in which explicit 
solutions for dynamical models on SF networks with in- 
finite size but finite cutoff were obtained [l3j]. In this 
case, an explicit dependence on k c was found, leading to 
a radically different behavior from that obtained for truly 
infinite SF networks, which then do not correspond to the 
analytical continuation to infinite k c . The true thermo- 
dynamic limit of SF networks corresponds to the double 
limit A^ — > 00 and k c — ► co, since keeping k c fixed leads 
to a finite degree second moment even for 7 < 3, and to 
a non SF network. The natural way to take this double 
limit is to allow k c to diverge with A~ but not faster than 
N 1 / 2 if the network is to be uncorrelated Hence 
any choice k c (N) ~ N 1 ^ with 2 < u> < co is legitimate 
and the way to reach the thermodynamic limit is not 
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unique. Therefore, in systems with anomalous finite size 
effects, depending on the way k c scales with system size, 
the final behavior as a function of N is modified, and will 
depend on the particular value of uj chosen. The hetero- 
geneous mean-field (MF) theory for CP turns out to dis- 
play anomalous scaling, therefore our findings clarify the 
discrepancies between MF theory and random-neighbor 
simulations recently reported in the literature [lfj] . More 
generally, our results indicate that the FSS analysis with 
the standard form (i.e. as a function of N alone) used in 
several recent works for SF quenched networks should be 
reconsidered. 

We start our case by considering the contact pro- 
cess on heterogeneous networks, which is defined as 
follows 3] ■ An initial fraction po of vertices is randomly 
chosen and occupied by a particle. Dynamics evolves 
in continuous time by the following stochastic processes: 
Particles in vertices of degree k create offsprings into their 
nearest neighbors at rate A/ k, independently of the de- 
gree k' of the nearest neighbors. At the same time, par- 
ticles disappear at a unit rate [Uj]. This model under- 
goes a continuous transition, located at a critical point 
A c , separating an absorbing phase from an active one 
with everlasting activity [15j|. The critical properties of 
CP on uncorrelated SF networks have been studied by 
means of heterogenous MF theory [7fl and a FSS the- 
ory Q. Heterogeneous MF theory predicts in the ther- 
modynamic limit a critical point A c = 1, independent 
of the network topology, a stationary particle density in 
the active phase given by p ~ A' 5 , with A = A — A c 
and (3 — 1/(7 — 2), and a density decay at criticality 
p(t) ~ t~ e , with 6 = 1/(7 - 2). On the other hand, 
Ref. [H assumes a standard FSS, depending only on N, 
which yields the particle density in surviving runs at crit- 
icality 
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for any k c > TV 1 / 7 . Apart from the unsolved question of 
the validity of heterogeneous MF theory for CP on un- 
correlated SF networks @, 0, EU , there is surprising evi- 
dence of a disagreement between the MF FSS exponents 
and simulations on a random neighbor version of SF net- 
works [lfj. In order to gain an understanding of finite 
size effects in the CP, we propose to focus on spreading 
experiments (lBj . i.e. in simulations starting with a single 
active site (p(t = 0) = po = 1/N) in which activity is fol- 
lowed until the systems decays to the absorbing state. In 
such experiments it is customary to measure the survival 
probability P(t), defined as the probability that activity 
survives up to time t. At the critical point this quantity 
scales as [15j) 



P(t) = t- 5 .f(t/t c ), 



(2) 



where the scaling function f(x) is constant for small val- 
ues of the argument and cutoff exponentially for x ^> 1, 



and t c is a characteristic time, depending usually on the 
size of the system. Standard homogeneous MF FSS the- 
ory predicts that at criticality the temporal cutoff scales 
as t c r~ N 1 / 2 [HI. To gain insight into the spreading ex- 
periment, we can map the dynamics at the MF level onto 
an effective one-dimensional diffusion problem. When a 
spreading experiment for the CP is performed the num- 
ber of active sites n(t) starts at 1 and at each time step 
can grow by 1, decrease by 1 or remain constant. Defin- 
ing the density of active sites pit) — n(t)/N, the diffusion 
process is defined in uncorrelated networks at criticality 
by the transition rates 
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(3) 

where the last term in Eq. ([3]) represents the probability 
that a random neighbor is empty, and Pkit) is the density 
of active sites restricted to nodes of degree k. Clearly 
pit) = /9fc(i)P(fc). The equation of motion for the pk 
at criticality is [7( 

d t p k it) = -p k it) + — [1 - Pk (t)]p(t). (4) 
(fc) 

Since all pk are expected to change in time slower than 
exponentially it is correct to perform a quasi-static ap- 
proximation [3] and set the l.h.s. of Eq. (g]) to zero, so 
that 



Pkit) 



kpjt)/(k) 
l + kp(t)/(k) 



(5) 



at sufficiently large times. Provided pit) >C (k)/k, which 
is always true in spreading experiments, the /c-dependent 
term in the denominator can be neglected, i.e. Pkit) — 
kpit)/ (k). This relation is very well verified numerically. 
Inserting it into the transition rates yields 



p^p- 



Pit), 



p^p^ 



pit) [1 - gp(t)] , (6) 



where g = (k 2 )/(k) 2 . Hence the density pit) of active 
states performs a biased, one-dimensional random walk, 
and the distribution Eq. ^ is the solution of the first 
passage time problem in p = 0, starting from the initial 
condition po = 1/N. Using the transitions rates ([6]) it 
is possible to measure the time t c with arbitrary preci- 
sion in trivial simulations of the random walk process. It 
turns out that this characteristic time depends explicitly 
on both the network size N and the network cut-off k c , 
through the form 
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(7) 



For 7 > 3, (non SF networks), g is a constant and Eq. ([7|) 
coincides with the usual MF result. For 7 < 3, on the 
other hand, g diverges with the upper cutoff k c of the 
degree distribution as fc^~ 7 , giving an anomalous form of 
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Since t c > t* for 7 < 3 and oj > 2, one has at the critical 
point 



t/(N/g) 



FIG. 1: Scaling plot of the survival probability of the CP on 
random neighbors SF networks with 7 = 2.5 and k c < N 1 ^ 2 . 
Minimum deg 2. Numerically we obtain 5 = 1. 



FSS, which separately and explicitly involves the system 
size N and the cutoff k r . i.e. 



(8) 



We have verified numerically, by performing spreading 
experiments on a random neighbors version of the un- 
correlated configuration model [IT]] with k c < N 1 ' 2 , that 
the surviving probability of the CP at criticality is well 
described by Eq. @ with a characteristic time t c given 
by Eq. ([7]), see Fig. Q] As we have discussed above, the 
cutoff k c is not fixed, and its divergence with network 
size can be chosen with a degree of arbitrariness. Impos- 
ing therefore k e (N) ~ N 1 ^ and expressing t c only as a 
function of N, this arbitrariness results in a dependence 
on the exponent ui, 



t ~ tv [1+(7 ~ 3)/ " i/2 . 



(9) 



This anomalous dependence of t c on both N and k c trans- 
lates into a density of active sites in surviving runs p s (the 
usual observable for FSS in non-equilibrium absorbing- 
state phase transitions) depending also separately on N 
and k c . To determine its detailed expression, it is cru- 
cial to realize that the solution of the MF Eq. ([4| for 
p(t) at criticality shows a crossover at a temporal scale 
t* ~ k~ c '~ 2 , in full analogy with the crossover occurring 
for diffusion-annihilation processes in SF networks [III]. 
The origin of the crossover can be traced back to Eq. |(SJ| . 
When there is a finite initial density po of active sites, the 
dependent term in the denominator of Eq. (|5|) cannot 
be neglected and hence, for t x < t < t* , the density de- 
cays as if the network size were effectively infinite, so that 
p(t) ~ i -1 /(T- 2 ) Q i x i s a microscopic time scale, inde- 
pendent of N and k c . As the density decreases, for times 
larger than t* , such that k c p(t*) = (k), the denominator 
of Eq. ([5|) does not effectively depend on k anymore so 
that for t* < t < t c , p{t) ~ (gt)^ 1 for any 7, yielding 
an exponent 9 = 1. The asymptotic value in surviving 
runs p s is reached for times larger than t c , so p s ~ p(t c )- 



gt c 
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Eq. (fT0|) shows that not only t c but also other observ- 
ables acquire an anomalous k c — dependence on finite SF 
networks, which translates into values of the exponent ra- 
tio (3/v depending on the particular u> considered. This 
result accounts for the numerical difference in the (3/v 
value found in random neighbor simulations with differ- 
ent values of uj [Iof . which can be better fitted [3] by 
means of Eq. (fTU]) , whose decay exponent is smaller than 
the one proposed in Refs. namely Eq. ([T]). It is 

also worth stressing that, since t c > t* , the true asymp- 
totic MF value for the temporal exponent 9 describing 
the decay of density is 1 for any 7 and any N < 00, and 
not the exponent 1/(7 — 2) predicted for 7 < 3 by the 
heterogeneous MF treatment for strictly infinite system 
size BB. 

This finite size separate dependence of scaling time at 
criticality on both the system size and the upper cut- 
off is not in fact a peculiarity of the CP in MF, but is 
actually found also in other processes taking place on 
SF networks. Sood and Redner [l8| have recently com- 
puted the time needed by the voter model dynamics tak- 
ing place on heterogeneous networks to reach a full con- 
sensus (ordered) state, which is given by r ~ N/g, and 
hence depends on both N and k c for 7 < 3. On the 
other hand, the distinct dependence on and k c is not 
the only possibility: There are systems for which the size 
scaling on heterogeneous networks is not anomalous and 
depends exclusively on N . Two examples are provided 
by the "link-update" [l9[ and "reverse" [2(| versions of 
the voter model. In the first case the time to reach full 
consensus is proportional to N [201 ] , with no dependence 
on network features whatsoever. For the reverse model 
instead such time is proportional to N(k)(k~ 1 ) [2(|. For 
all 7 > 2 these moments of the degree distribution are 
well-behaved and no dependence on k c arises. In this 
same spirit, it is possible to devise stochastic systems 
with continuous phase transitions that show standard 
finite size effects, depending only on the network size. 
Consider, for example, the variation of the CP proposed 
in Ref. [2l| (see also in which particles in vertices of 
degree k create offsprings in nearest neighbors of degree 
k' at rate A(/c£/)~ M , while particles still disappear at unit 
rate. For fi = 1, the heterogeneous MF theory for this 
model takes, in uncorrelated networks, the form 



d t pk(t) = -pk(t) 



Ap(t) 
(k) 



[1-P*(*)], 



(11) 



which has a threshold A c = (fc), and a stationary particle 
density in the active phase p ~ A — A c . It is easy to see 
that the random walk mapping of spreading experiments 
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7 > 3, and modified versions of the CP. Other models, on 
the other hand, show an additional explicit dependence 
on the upper cutoff of the degree distribution. In this 
last case the thermodynamic limit of an infinitely large 
network can be reached in different ways, depending on 
how k c is chosen to scale with TV, and different possible 
routes to the thermodynamic are thus possible. Under- 
standing what physical ingredients select which class a 
model belongs to is a challenging task for future work. 
R.P.-S. acknowledges financial support from the Span- 
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FIG. 2: Scaling plot of the survival probability of the modified 
CP in Ref. [21[ with fi = 1, on random neighbors SF networks 

with 7 = 2.5 and k c < N 1 ^ 2 . Minimum degree fc m i n = 2. — 
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FIG. 3: Scaling plot of the survival probability of the CP 
on random neighbors SF networks with 7 = 2.5 and k c = 
7V 1 /(7-i) Minimum degree fc min = 2. Notice that, for u> < 2, 
g displays very large sample to sample fluctuations. 



at criticality in this modified CP model is defined by the 
transition rates in Eq. (|6|) with g = 1, which correspond 
to a characteristic time scaling as t c ~ TV 1 / 2 , completely 
independent of the network cutoff k c . This result is fully 
confirmed by numerical simulations, see Fig. [2j 

Finally, let us discuss what happens in the standard 
CP when the cutoff scales as k c ~ N 1 ^ but with u < 2. 
This case includes the so-called natural cutoff ui = 7 — I 
(for 7 < 3) that occurs in the configuration model if one 
does not impose any constraint on the maximum degree 
of the vertices 
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In such a case, the presence of 
very large hubs invalidates the derivation of the Eq. (7). 
Figure [3] shows that this is in fact the case for 7 < 3 
and u> — 7 — 1 < 2, and it is then the proof that for 
random neighbor networks with k c > N 1 / 2 a different 
scaling form of t c must be considered. 

In conclusion, the evidence presented in this work leads 
to the conclusion that processes taking place on SF net- 
works may be divided in two classes. In the first class, 
finite size effects and FSS depend only on the number of 
nodes N. This class encompasses "reverse" and "link- 
update" voter dynamics as well as models showing con- 
tinuous phase transitions, such as the Ising model for 
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